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ABSTRACT 

We have found a new mean field solution in the BCS theory of superconductivity. This uncon- 
ventional solution indicates the existence of superconducting phase transitions of third order 
in thin films, or in bulk matter with a layered structure. The critical temperature increases 
with decreasing thickness of the layer, and does not exhibit the isotope effect. The electronic 
specific heat is a continuous function of temperature with a discontinuity in its derivative. 
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1 Introduction 

The BCS theory of superconductivity |]1| has been very succesful in explaining properties of a 
large class of simple superconductors in terms of just two experimental parameters, namely the 
squared electron-phonon coupling constant g and the Debye frequency ujd- The BCS theory in 
its simplest form is based on the following so-called BCS reduced grand canonical Hamiltonian, 

H = X!?k«Lak<7 -TtYI 0(^^d- I ?k \)&{huD- I <^k' l)4,+a-k,-a-k',-ak',+ (1) 

k,CT ^ k,k' 

where the quantities k and k' are wave vector variables, the quantity a = ± is the spin 
projection, and the symbols Oko- and aj^.^ are electron field operators. Furthermore, 

a-^-. (2) 

where // is the chemical potential and m is the electron mass. Finally, V is the volume of the 
system. The sum over k and k' in Eq. (P is restricted by the conditions |^k| < hu^, and 
|(^k'| < fi-^D, as indicated by the appropriate stepf unctions in Eq. (|ID. 

The Hamiltonian (|l|) represents by itself a great simplification of the net interaction be- 
tween the electrons. Even so, in BCS theory one makes an additional technical simplification 
by adopting merely the mean field approximation in calculating physical properties of super- 
conductors 0. These simplifications notwithstanding, the BCS theory predicts a number of 
relations, some of which are even independent of the phenomenological parameters g and uc, 
which are surprisingly well obeyed by a large class of superconductors [Q . 

As is well known, before the BCS theory was developed, Ginzburg and Landau had pro- 
posed a theory of superconductivity, currently known as the Ginzburg-Landau theory [Q. The 
equations of this theory were derived from the phenomenological Landau theory of second- 
order phase transitions . In accordance with the Landau theory, in all cases of second order 
phase transitions, one has to define the so-called order parameters, which characterise ordered 
structures of macroscopic fields appearing spontaneously in systems below a critical tempera- 
ture Tc. Thus the order parameters equal zero for T > Tc and acquire no n- vanishing values 
for T < Tc. Later Gor'kov demonstrated [^] by using the Green function method, that the 



Ginzburg-Landau theory is a limting case of the BCS theory, provided that the order param- 
eters are values of a complex function v^(x), which is proportional to the anomalous Green 
function < '0_(x)'?/'+(x) > of the electron field operators ^o-(x), evaluated in the mean field 
approximation. The expression for the function (/?(x) is thus the following, 

C 
V9(x) = C < V'-(x)V^+(x) >mf= ttYI! < ak-ak',+ >m/ expz(k + k') -x, (3) 

^ k,k' 

where < ... >„/ denotes the statistical grand canonical ensemble average with the Hamiltonian 
(|l]) in the mean field approximation. From Eq. (^) one sees that ordered superconducting struc- 
tures are characterised by the mean field correlation functions < ak,-ttk',+ >m/ corresponding 
to the Hamiltonian (|I|). Thus, the order parameters will be chosen to be related to these mean 
field correlation functions. 

A very important theorem concerning mean field solutions to the BCS theory of supercon- 
ductivity has been proved by Bogoliubov 0]. This theorem states that any mean field solution 
of the BCS Hamiltonian (|l|) becomes an exact solution in the thermodynamic limit. This 
means that the effects of quantal fluctuations about a given mean field give no contributions 
to thermodynamical potentials. This was also proven explicitly in [|| by a direct evaluation of 
the grandcanonical partition function corresponding to the Hamiltonian ([^). 

Although the BCS theory has existed for almost fourty years, apparently only one of its 
mean field solutions corresponding to the Hamiltonian (|l|) has been considered so far. To the 
best of our knowledge, there have been no attempts to explore the possibility of finding more 
than one mean field solution with the Hamiltonian (|I]). 

In this paper we report on our observation that there is at least one additional mean field 
solution to the system described by the Hamiltonian ([^). We describe this solution explicitly 
below. The superconducting state corresponding to the new mean field solution can appear at 
non-zero temperature only in quasi-t wo dimensional systems such as thin films or layered struc- 
tures of bulk material. Such layered structures are actually present in high Tc superconductors 
[0 . The new mean field solution described in this paper corresponds to a phase transition to a 
superconducting state with some rather unexpected properties. Specifically, it is a third-order 
phase transition which does not exhibit the isotope effect which is typical for the usual BCS 



theory. Thus, this phase transition cannot be described by the ordinary phenomenological 
Ginzburg-Landau theory, which describes phase transitions of second order. 

2 The mean field approximation and the gap equation 

In order to analyse the novel mean field solution to the BCS theory with the Hamiltonian (|lD 
we first express the interaction Hamiltonian Hi in (||) as a sum of two terms, 

Hj = H'j + H'l (4) 

where 



iij — 
and 



H'l = -— J2 0(^^D- I ^k \)<dihujD- I ^k' |)<5|ki,ik'| <^+aLk _a_k' -ak',+ (5) 

^ k,k' 



Hj = -— J2 ©(^^D- I ?k \)QihujD- I ^k' |)(1 - (^|k|,|k'|)ak,+alk -O-k' -ak',+ (6) 

^ k,k' 

In the interaction term defined by Eq. (^) we introduce new summation indices p and q 
related to k and k' as follows, 



k + k' k-k' 

P = ^-,q=^- (7) 

The innocent looking change of summation variables introduced in (|^) is very important, as it 
enables one to satisfy the condition | k | = | k' | in the sum in Eq. (|^) by any two perpendicular 
vectors p and q. This will be seen to be quite essential for the definition of order parameters 
below. 

The Hamiltonian (|lD can now be written as follows, 

H = ^Ckflkc.'^ko- - y^ <d{huJD- I ^p+q |)(5|p+q|,|p_q| ap+q^ + a_p_q _a_p+q _ap_q,+ + Hj (8) 

k,(T p,q 

Let us note that it is enough to retain only one step function in the equation (|) above, since 
^P+q = ^p-q in the sum in Eq. (|). 



We next introduce a general set of complex order parameters 6q,p and 6q p by the following 
definitions, 

Kp = 9 < a-(p+q) -«(p+q),+ > ^]p+q|,|p-q| , ^q.p = 9 < «(p+q), + al(p+q) _ > (^Ip+q] ,|p-q| (9) 

The order parameters 6q_p and 6* are thus enumerated by two orthogonal vectors p and q, 
which are restricted by the condition | ,^p+q | = | ^p_q |< hujD. 

We now define the following gap functions Ak and A^ in terms of the functions 6q^p and 



fe* p introduced above, 



Ak = :^ E ^,p5p-q,k, A* = i E 6;,p5p_q,k (10) 



p,q p,q 



Using the definitions (|]) and (|I0|) , we decompose the Hamiltonian (H) as a sum of two terms, 

H = H^f + Hfi, (11) 



where 



^™^ " 777 E ^q,P^-q-p + 

y p,q 



+ E {^k(ak,+«k,+ + 4 _ak ,-) - (A^a-k -ak,+ + A]^al^^a\_)j (12) 



k 

is the mean field Hamiltonian, and the term 



^n = -yJ2'^i^^D-^p+ci 



«{p+q), + "-(p+q) - < ^(p+q),+ ^-(p+q) - >. 



p,q 

(p-q) -a(p-q),+ - < a-{p-q) -a(p-q),+ > 



(^Ip+qMp-q I +Hj (13) 

describes quantal fluctuations about the mean fields 6q,p and 6* p, cf. Eq. (H). 

In view of the Bogoliubov theorem quoted above, the quantal fluctuation Hamiltonian (|13D 
has no macroscopic effects in the thermodynamical limit. By neglecting it in all expressions, 
i.e. replacing the Hamiltonian (§) with the mean field Hamiltonian (0) one gets all physical 
quantities in the so-called mean-field approximation. From now on we use this approximation 
throughout. 



The Hamiltonian Hmf, eq. (|T2|), can easily be diagonalized, and results then in the following 
partition function Zj^j, 

Zmf = rrexp[-/5iJ„/] 

= ^^P {-^ E ^^.P^-'^.p} {n {e-'Hl + e^^^)(l + e-^^--)) } , (14) 

where 

i?k = (e^ + |Akp)^ (15) 

is the quasiparticle energy spectrum. 

Evaluating the correlation functions in the mean field approximation, one then gets the 
following relations from the equations (|^), 

9 



^'^'P ~ 2Er. 



-Ap+q th (^/9i?p+q) • (16) 



^p+q 

The sets of equations (|TB|) for the determination of the functions 6q p and 6* p are closedly 
related to the so-called gap equations in the standard BCS-theory of supeconductivity, as will 
be seen presently. 



3 The solutions 6a d and b* ^ 

Vl)f q,p 

Using now equations (0) and (p!6D, one gets the following equation for the gap function Ak, 

^^=Y7^' ^^ ^h (^/5^P+q) Vq,k^|p-qMp+q|e(/la;^- | ^p+q I) (l^) 

V p,q ^^p+q V^ / 

We now revert to the summation indices k' = p + q, k" = p — q and carry out the summation 
over k" with the result 

^k = tIt E ^ th^E^' '^Ik'Ukl QifiLOo- I ek' I) (18) 

^1/ k' ^k' ^ 

The assumption that the gap function Ak' is a function of the magnitude | k' | only, simplifies 
the relation (pTsD substantially, as follows, 

^iViAk,,/3 



where 



iVi = EVMk|e(n^z5-ia'l) (20) 

k' 



The quantity A^i defined in Eq. (EDf) is the number of wave vectors k' with given length 



I k' 1 = 1 k |. The vectors k' have the components given below, 

, , (111 2-n 2n \ , , 

k = — rii,— n2,— na (21) 

where ni, n2 and 77,3 are integers, Li, L2 and L3 are the lengths of the edges of the box into 

which the system is enclosed. The condition | k' | = | k | defines the surface of an ellipsoid with 

semiaxes ai = ■^, 02 = ■^, and 03 = ■^, i.e. 

222 
ni; ni ni , , 

4 + ^ + ^ = 1 (22) 

Cfci dn ^"^ 

Thus the number of states A^i is in fact equal to the surface area of the ellipsoid (p2D which 
is determined by an elliptic integral. However for a very thin film, or layer with Li ^ L2, 
Li <^ L3 the number of states can be approximated by the surface area of a very thin disc of 
elliptical shape, with semiaxes 02 = -^ and 03 = ^^. Thus the number of states A'^i is given 
by the following expression, 

N,=27re^^ = ^eL2L, (23) 



Since the vectors k' in the sum ( |18D are restricted by the inequalities below, 

Sf - hiOD < <Sf + hujD (24) 

2m 

and since huD <^ £p we can replace the quantity k"^ in ( pSj ) by kp, i.e. 

iVi = ^k^L^Ls (25) 

ZTT 

By using the expression for the density of states of electrons with a single spin projection at 



the Fermi level A^o, 



we can express the ratio |^ as follows. 



2mkp , , 

^° - (2^ ^2^) 



^^ -N ^ ^ "^ (97) 



Using the last relation in Eq. (|T9|) we obtain the formula, 



Ak = ^^.gNoJ SpSi-^ th-(3E^ , (28) 

where 

is the lower bound of the energy spectrum of free electrons in a thin box of volume V = L1L2L3 



with L2 ^ Li and L3 ^ Li. From eq. (|28D follows that the energy spectrum E]^ of the 



quasiparticles in the superconducting state with A^ 7^ must satisfy the following relation, 

1 , 1 



1 = gNoJepSi— th-(3E^ , (30) 

which equation can be satisfied only if the temperature T is bounded from above by the following 
critical temperature Tc, 



Tc = -—gNo^SpSi . (31) 

ZKb 



The critical temperature Tc defined by eq. (31) does not depend on the Debye frequency 
Uo and hence it does not exhibit the isotope effect typical for the standard solution in BCS 
theory. The critical temperature also depends on the lowest bound Si of the energy spectrum 
of the electrons in the superconducting material. For bulk materials £^1 ^ so this kind of 
superconductivity does not exist in bulk samples. However for thin films the lower bound Si is 
nonzero and given by Eq. (^). The critical temperature increases with decreasing thickness 
of the layer. 

The solution to eq. ( ^OD for the energy spectrum E'k in the range of k for which Ak 7^ 0, is 
a quantity ti{T) say, which is dependent on the temperature T only, 

i?k = (C^+|AkP)^ = r^(r) (32) 

One cannot obtain an explicit expression for the solution ri{T) = E'k of eq. (BO) for general 



values of T in the range T < Tc but has to resort to numerical methods for these general values 
of T. However in the limiting cases T <^Tc and Tc — T<ti Tc, respectively, one readily obtains 
the following formulae from eq. (|5D|), 



77(T) ^770(1- 26-'^"") ;^<rc, (33) 



and 



where 



The ratio 



rj{T)^V3r]o{l-^y ■,T,-T<t:T,, (34) 



Vo = V{T = 0) = 2kBT,. (35) 



''^ 2 (36) 



ksTc 

is a universal constant, independent of the physical properties of the thin layer under consid- 
eration. 



From eq. (|32D we get the following expression for the gap function Ak, 

Ak=(^^-e9'e(r^-|a|). (37) 

Combining eq. ( pTj ) with (|T6|) one finally obtains an explicit expression for the order parameter 



"q.P' 



\,. = ^-.,P = 4^Jv' 



'^^^'^-'^''^ 



0(^-||^(q' + p')-^Fl)- (3^ 



The last two formulae, (^) and (^), are of course consistent with the definition (p!oD. The 
formulae (|30| ) - (^) are valid under the following condition, 

ri{T) < huD (39) 

because for |^k| > ^^d the quasiparticle energy spectrum i?k must reduce to the energy spec- 
trum of free electrons. Combining the relation (|35D with the inequality (p9D one obtains. 



r.<^^5T. (40) 

where T^ is the Debye temperature. For the case under consideration we thus have an upper 
bound for the critical temperature. Critical temperatures T^ in the vicinity of this upper bound 
have in fact been observed in layered structures of high Tc superconductors . 

The expressions (^7|) and (^) indicate unusual properties of the corresponding supercon- 
ducting state. Namely, not all of the order parameters bq^p with p • q = acquire spontaneously 

8 



non-vanishing values at the same time when T < Tc. For T = Tc — only one order parameter 
is non-vanishing at the opening of the energy gap r]{T = Tc — 0) = 0+. Decreasing the temper- 
ature from the value T^ the gap in the energy spectrum r]{T) opens up more and more, whence 
more and more order parameters acquire non- vanishing values for T kT^. 

4 Thermodynamical properties 

We next calculate thermodynamical properties related to the novel solution analysed above. 
Using the standard method we express the grand canonical potential VLgiT) of the super- 
conducting state by the formula 

2, 



6 Jo 



d 
drj' 



1]' 2 



(41) 



where fi„(T) is the normal-state grand canonical potential of the system if it were in a normal 
state at temperature T. The integral (^l|) requires numerical analysis for general values of T in 



the appropriate range. However, similarly as in the previous cases (0) and (p^) the following 
limiting cases are readily obtained. 



and 



Qs{T)-nn{T) = -^VNovl 



4VS 



4 \tJ 



,T<t:Tc 



n^iT) - n^{T) = -^vNoT]^, (^ ~ T ) ' ' ^<= - ^ ^ ^^ 



(42) 



(43) 



Using the formulae ( ^2]) and ( ^31) one can calculate the electronic specific heat CsiT) and the 
critical magnetic field Hc{T) with the following results. 



C,{T) = AVNoVoksiPvofe-^'^' .T <^% , 

2n 



57r 



Hc{T) = [—N,] r^o 



Likewise, in the limiting case corresponding to eq. ( p4D one obtains. 



TT^ fTy 



(44) 
(45) 



^^(^)-^"(^)-6^^1-^V,T.-T«T. 



CniT) 



TT^ 



Tr 



(46) 



as well as 



1 5 



H.iT) = f^iVo j r^o (l - ^) ', T. - T « T, (47) 



where 

27r2 



Cn{T) = ^VNoklT (48) 



is the normal-state electronic specific heat. 



In the low-temperature region T <^ Tc the formulae (^) - ( ^Tf ) do not yield results which 
differ appreciably from the standard solution 0. However, near the critical temperature, i.e. 
for Tc — T <^ Tc the differences are crucial. The formulae (^), (^61) and (^) yield critical 
exponents which differ from those known for the standard solution. The electronic specific 
heat given by eq. (^6D is a continuos function of T at T = Tc, however with a discontinuouss 
derivative at T = Tc, in contradistinction to the standard case. The novel mean field solution 
considered here gives rise to a superconducting phase transition of third order, and can therefore 
not be described by the ordinary macroscopic Ginzburg-Landau theory of superconductivity, 
which is relevant for the phase transitions of second order. 

5 Summary 

The unconventional superconducting state which we have found within the mean field treatment 
of BCS theory can appear in thin films, or in layered structures in bulk material. Therefore, the 
existence of this novel feature of BCS theory can have relevance for high T^ superconductors or 
for the explanation of the fact that even some simple metals can become superconducting as 
thin films, but not as bulk material. The fact that the BCS Hamiltonian contains the possibility 
of two different superconducting structures, namely an unconventional superconducting state at 
a critical temperature Td, and the well-known superconducting state at a critical temperature 
Tc, offers the possibility of observing consecutive superconducting phase transitions in thin 
layered structures. 
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